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Abstract In view of the recently observed discrepancy of theory and experiment for muonic hy- 
drogen [R. Pohl et ai, Nature vol. 466, p. 213 (2010)], we reexamine the theory on which the 
quantum electrodynamic (QED) predictions are based. In particular, we update the theory of the 
2P-2S Lamb shift, by calculating the self-energy of the bound muon in the full Coulomb+vacuum 
■ polarization (Uehling) potential. We also investigate the relativistic two-body corrections to the 

vacuum polarization shift, and we analyze the influence of the shape of the nuclear charge dis- 
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tribution on the proton radius determination. The uncertainty associated with the third Zemach 
moment (f 3 ) 2 in the determination of the proton radius from the measurement is estimated. An 
updated theoretical prediction for the 2S-2P transition is given. 
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Section 1: Introduction 



1 Introduction 

The high accuracy of quantum electrodynamic (QED) predictions and the precise spectroscopy of simple 
atomic systems allow for the determination of fundamental constants like the Rydberg constant i?^ from 
the hydrogen spectrum, a from the helium fine structure, and the electron mass m e from the g factor of 
hydrogenlike ions. In all these cases, nuclear structure effects arc small or can be eliminated. On the other 
hand, from a comparison of the theoretical isotope shift of optical transitions with experimental values for 
such atoms as He, Li, and Bc + , the nuclear charge radii of the heavy, unstable isotopes, in comparison to 
the stable ones, have been determined with high accuracy [TJ[2]- The most elaborate measurement [3] of 
the hydrogen-deuterium (H-D) isotope shift of 1S-2S transition gives a deuterium-proton mean square 
charge radius difference accurate to 2 parts in 10 5 . In some cases, the excitation of nuclei by the orbiting 
electron, the so called nuclear polarizability correction, is a significant effect and has to be included in 
the theoretical predictions. For example, the 1S-2S transition in D is affected by about 20 kHz due to 
the deuteron polarizability [3] , while the current experimental precision for the H-D isotope shift is below 
100 Hz (see Ref. 0). 

Bound muons penetrate the electron charge cloud and undergo transitions close to the atomic nucleus, 
partially screened by the remaining electrons. The determination of absolute charge radii from the muonic 
transitions is an established technique for the determination of nuclear radii [SJE|. Also, the CODATA 
values [7] for the proton and deuteron nuclear radii are mainly determined by an analysis of the 23 most 
accurately measured transitions in these systems [5J. Both atomic hydrogen as well as bound muonic 
systems are essentially two-body bound systems, with a comparatively light orbiting particle (electron 
or muon) and a heavier nucleus. Still, the combined evaluation of relativistic, QED and two-body effects 
remains difficult, and the presently achieved accuracy is the result of decade-long work of a number of 
physicists. Since the original calculations of Lamb shift by Bethe, later by Feynman and others, theory 
has been worked out to a very high precision. 

For hydrogen, present limitations come from the inaccuracy in the two-loop electron self-energy and 
amount to about 1 kHz for the 15 Lamb shift. This constitutes a relatively small uncertainty in comparison 
to the 1 MHz shift due to the proton finite size r p and allows for a determination of r p with an accuracy 
of 8 parts per thousand, namely 

r p = 0.8768(69) fm (Refs. 0H]) . (1.1) 

This value is in excellent agreement with the recently obtained proton radius from electron scattering [5] , 
which yields a value of r p = 0.879(8) fm. Here, we have added the statistical and systematic uncertainties 
given in Ref. [9] quadratically. The muonic hydrogen value 

r p = 0.84184(67) fm (Refs. [TO]) (1.2) 

is in significant (5.0 a) disagreement with the CODATA and (4.6 er) with the electron scattering value. 
This proton radius significantly disagrees with a number of investigations, including determinations of the 
proton radius based on previous scattering experiments ("world scattering data") reported in Refs. [111 - 
H3]. 

For atomic hydrogen and more generally for electronic bound systems, the energy shift due to the finite 
nuclear size is proportional to the mean square charge radius (r 2 ), and the contribution from higher 
moments such as the convoluted, third-order Zemach moment (?" 3 ) 2 is negligibly small [2], only about 
50 Hz, and similarly, proton polarizability effects amount to only about 100 Hz (see Refs. [T5"Hl7| ). Both 
of the latter effects thus have only a negligible effect on the proton radius determination from electronic 
hydrogen spectroscopy (we sometimes refer to atomic hydrogen as "electronic hydrogen" in the current 
work in order to uniquely distinguish the system from muonic hydrogen). 

Indeed, the situation is different in muonic hydrogen. The muon Bohr radius is about 200 times smaller, 
and the energy levels are significantly affected by the proton finite size. The contribution from higher 
order moments of the nuclear charge distribution and from the proton polarizability, although still small, 
have to be accurately estimated in order to obtain a reliable proton charge radius. If one refrains from 
determining a new proton radius from the recent measurement [TOJ of the 2S(F=1) 2P 3 / 2 (F = 2) 
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transition and uses the CODATA proton radius in order to predict the transition energy, then one finds 
that the experimental result reported in Ref. [10j deviates from QED theory by about 0.31 meV. 

A reinvestigation of the theory of the 2S-2P Lamb shift transition in muonic hydrogen is thus indicated, 
especially because the theory of muonic atoms is not free of surprises [J8]. We proceed as follows. In 
Sec. [U we first discuss the size of the observed discrepancy in relation to the relativistic and QED 
effects which describe the energy levels of muonic hydrogen. We then continue to verify the current 
status of theory on the basis of the two-body Breit Hamiltonian and evaluate relativistic effects and 
quantum electrodynamic corrections (vacuum polarization and self energy, and combined effects) which 
are relevant on the level of the current disagreement. We then continue to calculate a few hitherto 
neglected corrections to theory (Sec. [3]). Special attention is devoted to the calculation of the Bcthc 
logarithm in the full Coulomb+Uehling potential, which amounts to a nonperturbative treatment of the 
vacuum polarization correction to the self energy. A final update of theory and a rccvaluation of the 
proton radius is presented in Sec. |4j Conclusions are reserved for Sec. [5] 



2 Verification of Theory 

2.1 Size of the Discrepancy 

In view the recently observed discrepancy of theory and experiment in muonic hydrogen [10) , we attempt 
to verify the theory [T9"H2"7] used for the theoretical evaluation of the recent measurement [TU] . We base 
our considerations on the extensive list of corrections presented in Tables 1 and 2 of the supplementary 
information included with the recent paper [10], whose entries in turn are based on previous theoretical 
work [T5H27j. 

The authors of Ref. [TU] report on a measurement of the energy interval 

AE = E (2P 3/2 (F = 2)) - E (2S 1/2 (F = 1)) (2.1) 

in muonic hydrogen. Summing up all known effects, the authors of Ref. [lOj use the following theoretical 
prediction 

E t h = ^209.9779(49) - 5.2262^ + 0.0347^^ meV. (2.2) 
The experimental result is 

P exp = 206.2949(32) meV . (2.3) 
One thus infers the root- mean-square proton radius given in Eq. (|1.2[) . 

On the other hand, one can also interpret the result of the measurement [TU] as a disagreement of 
theory and experiment in an important quantum electrodynamic measurement. Namely, if one uses the 
CODATA recommended proton radius given in Eq. (|1.1[) and evaluates theory according to Eq. (|2.2p . 
then one obtains a theoretical prediction of E t h = 205.984(63) meV which deviates from the experimental 
result given in Eq. (|2.3[) by 5.0 a. As already mentioned, an additional, hypothetical effect that would 
shift theory by (roughly) +0.31 meV might thus explain the discrepancy of the proton radius given in 
Eq. |L2]) and the 2006 CODATA value of r p = 0.8768(69) fm. 

The 2S\/2(F = 1) 2P 3 / 2 (F = 2) transition frequency is the sum of three contributions: 

• (i) the 2S-2P 1/2 Lamb shift, 

• (ii) the 2P 1 / 2 -2P 3 / 2 fine-structure splitting, and 

• (iii) the 2S and 2P 3 / 2 hyperfine structure. 

On the level of the theoretical-experimental discrepancy (0.31 meV), these contributions can be broken 
down into a set of well-established, conceptually simple relativistic and QED corrections, as detailed 
below. 
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2.2 Reduced Mass Dependence 

Our goal is to verify the theory of the 2S (F—l) <=?■ 2P 3 / 2 (F=2) transition energy in muonic hydrogen on 
the level of the reported discrepancy [10] , i.e., 0.31 meV in energy units. The unperturbed nonrelativistic 
Schrodinger Hamiltonian of the bound muon-proton system is 

^0 = /^--, (2.4) 

and the corresponding unperturbed binding energies of bound states in muonic hydrogen are given by 

a 2 m R 

E ° = ~ 2^T ' (2 ' 5) 

where tur = m M /(m^ + tojv) is the reduced mass of the system (we use units with h = c = eo = 1). 
Here, n is the principal quantum number. The proportionality to the reduced mass (not to the muon 
mass) follows from an elementary separation of the two-body Hamiltonian describing the proton and 
muon into relative and center-of-mass coordinates [28]. Relativistic corrections enter at relative order a 2 , 
i.e., at order a. A mR. Their reduced mass dependence is captured in the two-body Breit Hamiltonian |29j . 
For muonic hydrogen, the mass ratio of orbiting particle and nucleus is 

~m~ p ~ 8.880 243 ... ~ 9 ' ( ' ' 

and therefore not necessarily small compared to unity. The one-particle Dirac equation, which is valid 
only in the limit m^/m p —> 0, therefore cannot be used as a good approximation for the calculations of 
the relativistic effects. 

Vacuum polarization effects shift the muonic hydrogen spectrum at relative order a, i.e., at order a 3 mR. 
Although the vacuum polarization corrections are thus larger than the relativistic corrections, we start 
our discussion with the latter effects as they constitute the most natural extension of nonrelativistic 
atomic theory without held quantization. 



2.3 Relativistic Effects and Fine Structure 

The two-body Breit-Pauli Hamiltonian for the muon-proton system, including the reduced mass correc- 
tions but without the anomalous magnetic moment of the muon, reads as follows, 

p 4 p 4 a ,( 6 lj r % r j \ ■ 

8 8 TOp I nip \ r r° / 

+ ¥ (<•» + iii + dy rfw + ■ 9 "*m 

8 m p r 3 \ r 2 J 4 r 3 \ m 2 m M m p m 2 

The proton g factor is g p = 5.585 694 712(46) (see Ref. [7]). The anomalous magnetic moment of the 
proton comes from its internal structure. In the leading QED approximation, one assumes a point proton, 
minimally coupled to the electromagnetic field with its physical value of the magnetic moment which 
determines the hyperfinc splitting. Strictly speaking, this phcnomcnologically inserted intrinsic proton 
magnetic moment leads to a QED theory which is not renormalizable. The situation is different for the 
anomalous magnetic moment of the muon, which is described by vertex corrections within QED theory. 
Lacking a convenient way to describe ab initio the inner structure of the proton with its anomalous 
magnetic moment, the phenomenological insertion of g p into the proton-spin dependent terms in the 
Breit-Pauli Hamiltonian is the most obvious treatment in theoretical calculations adopted so far in the 
literature. 



j - fx p. (2.7) 
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As already mentioned above, the 2S 1 / 2 (P = 1) •<=> 2P 3 / 2 (P = 2) transition frequency is the sum of 
three contributions: (i) the 2S' 1 / 2 -2P 1 / 2 Lamb shift, (ii) the 2P]y 2 -2P 3 / 2 fine-structure splitting, and (hi) 
the 2S\/2 and 2P 3 / 2 hyperfine structure. The 2S\/2 and 2P 3 / 2 hyperfine structure splittings, as well as 
the 2P 1 / 2 -2i- > 3/2 fine-structure splitting can directly be determined from the Breit-Pauli Hamiltonian 
given in Eq. (|2.7p . On the level of the theoretical-experimental discrepancy (0.31 meV), we find that the 
neglected terms beyond the Breit-Pauli approximation arc already too small in magnitude to explain the 
discrepancy on the basis of either the fine-structure or hyperfine-structure intervals. 

We start with the 25 hyperfine structure, which is exclusively given by the terms proportional to a p ■ 
5 3 (r) in the Breit-Pauli Hamiltonian. It is sufficient to evaluate matrix elements of this operator, with 
nonrelativistic spinor wave functions that are solutions of the Schrodingcr Hamiltonian of the system 
(with reduced mass mn). The result is 

E*£(2S) = (2S(F=1)\H BP \2S(F = 1)) - (2S(F=0)\H BP \2S(F = 0)) 

w 9p a 4^k_ = 22.805 meV. (2.8) 

6 trip 

After a consideration of the proton structure (Zemach), and QED radiative corrections, one obtains (see 
Ref. [25]) a theoretical prediction of 

£g™ (2S) = 22.8148(78) meV. (2.9) 
For the 25(^=1) sublcvel, the difference 

- A (X Q ™(25) - Ejg(2S)) = 0.0023 meV (2.10) 

already is two orders of magnitude smaller than the experimental-theoretical discrepancy of 0.31 meV. 

The 2P]y2 hyperfine splitting (energy difference of the F = 1 and F — states) is not needed for the 
analysis of the experiment of Ref. [I0]- Nevertheless, it is interesting to note that from Eq. (|2.7[) , we have 

E%Z{2P 1/2 ) = (2P 1/2 (F = l)|# B p|2Pi /2 (F = l)) - (2P 1/2 (P = 0)|P B p|2P 1/2 (P = 0)) 
_ g p a 4 m\ ( g p -l 



.1 + ^ ^ = 7.953 meV. (2.11) 

18 m M Trip \ 2g v m p J 

Including the anomalous magnetic moment of the muon, one obtains a result of |19j 

£ hfe D ( 2 - P i/ 2 ) = 7.963 meV, (2.12) 

which differs from P^(2P!/ 2 ) only on the level of 0.010 meV. 

The 2P 3 / 2 hyperfine splitting (energy difference of the F = 2 and F = 1 states) is obtained from the 
Breit-Pauli Hamiltonian as 

?hfe ( 2P 3/2) = (2-T3/2 K r = 

g p a Ar m? R ( 5 g p — 1 m 



pBP(2P 3/2 ) = (2P 3/2 (P = 2)|P B p|2P 3/2 (P = 2)> - (2P 3/2 (P = l)| J ff B p|2P 3/2 (P = l)) 



,1 + 7— = 3.392mcV. (2.13) 

45 rrifj, m p \ 4 g p m p J 

The full QED result is minimally displaced by less than 0.001 meV (see Ref. [27]), 

^i?™ ( 2 ^3/2) = 3.3926 meV, (2.14) 

The Breit-Pauli Hamiltonian also leads to an off-diagonal coupling of the 2P]y 2 (P = l) and 2P 3 / 2 (P = 1) 
sublevels [19], which evaluates to 

^ BP (2P) = (2P 1/2 (P = l)|P B p|2P 3/2 (P = l)) 



/2 \r =i)\n B p\^-r 3 / 
9 P a 4 m 3 R ( g p - 1 m M 



144 V2 ffift uip \ g p m p 



-0.796 meV, 
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in excellent agreement with the value given in Eq. (85) of Ref. [T9"] . 

Finally, the energy difference of the hyperfine centroids (the fine-structure splitting) of the 2P states is 
readily evaluated as 



(2P) 



(2P 3/2 \H BP \2P 3/2 ) 



2P 1/2 \H BP \2P 1 



32 m M m p 



1 



= 8.329 meV. 



This result differs by only 0.023 meV from the full QED result of [TW2"2]. 

E QJ^(2P) = 8.352 meV. 

If we define the states 

|1) = |2P 1/2 (F = 0)), |2) = |2P 1/2 (F = 1)), |3) = |2P 3/2 (F = 1)) 



(2.15) 



|4) = |2P 3/2 (P = 2)) 



and the matrix elements 

A/2 = £ hf8 (2P 1/2 ) 



V(2P) , (3 3/2 = E hfs (2P 3/2 ) , f = E Ms (2P) , 



(2.16) 



(2.17) 



and the zero point of the energy scale to be the hyperfine centroid of the 2Pi/ 2 levels, then the Breit-Pauli 
Hamiltonian in the 2P state manifold assumes the following matrix form M B p- 



M 



V 



[A/2 









J 01/2 





l ! 

;A»/2 + / 








03/2 + / 



(2.18) 



The off-diagonal elements v lead to admixtures to the |2P 1 / 2 (P = 1)) levels from the |2P 3 / 2 (F = 1)) 
levels and vice versa, and to a repulsive interaction as for any coupled two-level system. In agreement 
with this general consideration, a diagonalization of Mbp immediately leads to the conclusion, that the 
|2Pi/ 2 (P=l)) is lowered in energy by 

A = 0.145meV, (2.19) 

whereas the |2P 3 / 2 (P = 1)) energy is increased by A. This is in full agreement with Ref. [19] . 

In the current derivation, we have not distinguished the magnetic projections of the hyperfine sublevels. If 
these are included, a 12-dimensional matrix is obtained from the singlet |1), the three magnetic sublevels 
of states 1 2) and |3), and the five magnetic sublevels of state |4). However, we have checked by an 
explicit calculation of the Hamiltonian matrix, using angular momentum algebra [3THl3"T] , that the energies 
obtained by diagonalizing the full Hamiltonian matrix are the same as those obtained from (|2.18p . and 
that, in particular, the state |2P 3 / 2 (P = 2)), which is so important for the experiment |10j . remains 
uncoupled from the other hyperfine and fine-structure levels. 



2.4 QED and Lamb Shift 

The theory of the 2P 1 / 2 -2«S' 1 / 2 Lamb shift in muonic hydrogen is surprisingly simple. One here speaks 
of the 2P-2S Lamb shift because the 2P level is energetically higher than 25, in contrast to elec- 
tronic hydrogen where the situation is opposite, and the 2S-2P Lamb shift is observed. The one-loop 
vacuum-polarization (Uchling) correction due to virtual electron-positron pairs gives the main contribu- 
tion, and already the first-order effect (in perturbation theory) accounts for 99.5 % of the Lamb shift, 
or 205.0074 meV. The second-order perturbation theory effect contributes 0.1509 meV, and the two-loop 
vacuum polarization gives a shift of 1.5081 meV. We note that in the 1970s, an error in the evaluation 
of the two-loop vacuum polarization [32] has led to a discrepancy of theory and experiment for heavy 
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muonic ions |33j . The discrepancy was reduced after the error was discovered |34H36j . This observation 
implies that a careful verification of the vacuum polarization effects appears worthwhile in the current 
situation. 

Finally, the one-loop self-energy of the bound muon decreases the 2P-2S Lamb shift and amounts to 
— 0.6677 meV (the sign is natural as the self-energy effect is the dominant effect in electronic hydrogen, 
where the sign of the entire Lamb shift is reversed compared to muonic hydrogen, due to this term). The 
four mentioned contributions to the Lamb shift, which are recalculated and verified below, already account 
for 205.9987 me V. The total QED result without the nuclear finite-size effect, according to Ref. [10], is 
206.0573 meV. The difference, 0.0586 mcV, is already much smaller than the observed discrepancy of 
0.31 meV, and there are no cancellations among conceivably large, further QED effects not accounted for 
in the above list of dominant vacuum-polarization and self-energy effects. This illustrates that the theory 
of the Lamb shift in muonic hydrogen rests, to a large extent, on a very well established subset of QED 
effects. 



2.5 One— Loop Vacuum Polarization 

The one- loop vacuum polarization effect is described by the Uehling Potential V vp (r) which was evaluated 
already in 1935 (Ref. [37]). Muonic hydrogen is different from electronic hydrogen because the classical 
orbit of the bound muon is much closer to the proton than for the electron. It is thus useful to analyze the 
asymptotic behavior of the Uehling potential for large and small distances from the nucleus. We define 
the scaled coordinate p = amur, where a is the fine-structure constant, and ttir is the reduced mass of 
the bound system. The Coulomb potential scales as 



The one-loop Uehling potential is 



V{r) 



a a m r 
r p 



(2.20) 



Kp(r) 



TXp 



due 



3u 4 



(2.21) 



where x 



0.73738368 . . . for muonic hydrogen. In atomic hydrogen, x 



^> 1, and the 



Uehling potential can be approximated by a Dirac S, but this is not the case for muonic atoms. For small 
distances (p—t 0), one finds, in agreement with Ref. [34) . 



K P (r) 



o?m r 

7T p 



\{Hp*)+1e)-\px + \ 



O(p) 



(2.22) 



As compared to the Coulomb potential, only a logarithmic divergence is added. For large distances, 
p — > oo, by contrast, the Uehling potential decreases exponentially, 



K P (r) 



a 3 m r c _ 2px 



1 



2!) 



2225 



4^5/23,3/2 64^7/23,5/2 2048p 9 / 2 a; 7 /2 



O 



1 



,11/2 



(2.23) 



which leads to rapidly convergent radial integrals. On the Bohr radius scale of muonic atoms as measured 
by the scaled coordinate p, the Uehling potential goes as ln(p)/p for p — > 0. This is more singular than the 
Coulomb potential, but only by a logarithm. The smooth behavior excludes conceivable nonperturbative 
effects which could be expected for a singular behavior near the origin. 

The question then is whether one should start the evaluation of the one-loop vacuum polarization cor- 
rection from the nonrelativistic Schrodinger wave function or from rclativistic Dirac theory. The latter 
approach is indicated for highly charged ions, where relativistic effects dominate over the reduced-mass 
corrections. Muonic hydrogen, on the other hand, is a light system, and the parameter characterizing the 
reduced-mass corrections (m^/nip) is much larger than the parameter characterizing the expansion into 
Coulomb vertices (Za = a, where Z = 1 is the charge number of the proton). We therefore start from 
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the Schrodinger Hamiltonian and add reduced-mass effects later on the basis of a radiatively corrected 
Breit-Pauli Hamiltonian. 

Up to now (Refs. [HOOH]), ^he vacuum-polarization corrections have been evaluated using an analytic 
approach, with analytic representations of the reduced Green functions (sec Ref. |19j for relevant for- 
mulas). Here, in order to check for any conceivable calculational errors in the numerically dominant 
vacuum polarization corrections, we choose a different approach. Namely, because the singularity of the 
Uchling potential at the origin is only logarithmic, the combined Coulomb+Uehling potential is eligible 
for a numerical solution of the Schrodinger equation on a lattice |38j . In this approach, one may choose 
to evaluate perturbative terms using unperturbed Schrodinger eigenstates (Green functions are repre- 
sented by the pseudospectrum of states obtained from the finite lattice) , or one may choose to perform 
a nonpcrturbativc solution of the Schrodinger equation, for the combined Coulomb+Uehling potential. 

Performing a matrix element calculation on the lattice composed of an exponential grid (sec Ref. [39] ) 
and observing the apparent convergence on lattices with more than N = 300 grid points, and increasing 
TV in steps of AN = 50, we confirm the results, 



5E {1) = (2P|V vp |2P) - (25|y vp |25) = 205.0074meV 



in first order and 



5E^ 2 



2P 



U 



vp 



(E 



Ho)' Vvp 



2P 



25 



(E 



Ho)' Vvp 



2S ) = 0.1509 me V 



(2.24) 



(2.25) 



in second order. Here Ho and Eq are the Schrodinger Hamiltonian of the system and the binding energy 
given in Eqs. (|2.4j) and (|2.5j) . and [1/(Eq — Ho)'] denotes the reduced Green function. 

The Hamiltonian of muonic hydrogen, including vacuum polarization, is 



H 



Ovp 



P- 



2 ITLfj 



+ K P (r) , 



(2.26) 



For the numerical calculation, we add the Uehling potential to the Schrodinger Hamiltonian and diag- 
onalizc the Hamiltonian Povp on a pseudospectrum of states that is generated by a finite, exponential 
grid [38l[39]. We use the exact reduced mass of the system, which cannot be done consistently when 
using the Dirac equation, and a variable coupling parameter x = a / n > f° r various values of x, in order 
to control the convergence of the calculation. We find a smooth dependence on x, which can be fitted 
to excellent accuracy by a power series in \- The first two terms are consistent with the results derived 
in Eqs. (|2.24[) and (|2.25[) . For \ = a/n, which is the physically relevant coupling parameter for muonic 
hydrogen, our result for the total 2P-2S electronic vacuum polarization contribution to the Lamb shift of 
205.1584(1) meV is in excellent agreement with the sum of the one-loop term (205.0074 meV), the iterated 
one-loop term (second order, 0.1509 meV) and the third-order perturbation theory effect of 0.00007 meV 
as discussed in Refs. |401|4T1. 



2.6 Two— Loop Vacuum Polarization 

The two-loop vacuum polarization correction corresponds to the diagrams in Fig. [TJ The expression for 
the diagrams was first derived by Kallcn and Sabry in 1955 (Ref. |42j). and the corrections therefore carry 
their name. The two-loop calculation leading to their derivation is non-trivial, and it is thus imperative 
to clarify the status of the two-loop potential in the literature. A clear exposition of the derivation is 
not only given in volume III of the monograph [43], but the two- loop effect was also recalculated and 
later used as input for the three-loop corrections to the electron g factor; explicit results are indicated in 
Ref. [44]. Finally, the result has later been generalized to non-Abelian gauge theories, see Refs. [45h47] . 

An evaluation of the two-loop potential with nonrelativistic Schrodinger wave functions, using the two- 
loop potential K/p given in Refs. [H]|35] then leads to the result 

SE KS = (2P|U v ( p 2) |2P) - (2S|K ( p 2) |2S) = 1.5081 meV. (2.27) 
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Figure 1: Two- loop vacuum-polarization diagrams. 



For reference, it is quite instructive to divide the calculation into an evaluation of the loop-after-loop 
correction in the leftmost diagram of Fig. [T] (which gives a contribution of 0.25 meV to the 2S-2P Lamb 
shift), and the "true" two- loop effects diagrammatically represented in three rightmost diagrams in Fig. [I] 
which contribute 1.25 meV. The total value of 1.50 meV for the two-loop effect (T9"ll2"2] is thus confirmed. 



2.7 Muon Self-Energy 

The self-energy shifts of the 25* and 2P 1 / 2 levels are given by 

This expression has been derived originally in the early days of quantum electrodynamics |49| . and the 
effect of relative order a has been derived in Refs. [SDIEI]. The coefficients agree with very precise 
numerical investigations of the one- loop self-energy |52j . in the domain of low nuclear charge numbers, 
and an extensive list of Bethe logarithms lnfco(nP) has been given in Ref. [53]. 

There is a further effect due to virtual muon-antimuon pairs that modify the vacuum polarization poten- 
tial, 



The reduced-mass dependence of the corrections listed in Eqs. (|2.28p and (|2.29[) has been analyzed in 
Ref. [51]. 

The self-energy effect on the 2P state is given by 

E SE (2P 1/2 ) = {-\^L - \ \nk (2S)\ . (2.30) 

The first term in curly brackets is due to the anomalous magnetic moment of the muon. Its reduced-mass 
dependence, which is different from that of all other terms, follows from the Breit-Pauli Hamiltonian if 
one corrects the muon-spin dependent terms by the muon anomalous magnetic moment |54j . Otherwise, 
the coefficients in Eq. (|2.30[) have been verified against very precise numerical calculations [55] of the 
one-loop self-energy. There is no additional vacuum-polarization effect for the 2Pi/ 2 state to the order 
under investigation (a 6 mfl). 

The final result for the self-energy and muonic- vacuum polarization shift is 

SE SE = E SE (2P 1/2 ) - E VP (2S 1/2 ) - E SE {2S 1/2 ) = -0.6677 meV. (2.31) 
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p p p p 

Figure 2: Electronic vacuum polarization corrections to the muon self- 
energy. 

Finally, the contributions given in Eqs. (|2T2"4"]) . (|2T2"5|) . (|2~2"7]) and ([2"3T]) add up to 

5E = SE {1) + 5E {2) + SE KS + 8E SE = 205.9987 mcV , (2.32) 

All further QED corrections listed in Table 1 of Ref. [TU] (supplementary material) add up to 0.0586 meV, 
which gives a total QED result of 206.0573 meV. Both individually as well as collectively, the remaining 
effects are thus too small to explain the observed discrepancy of 0.31 meV. 

3 Advancing Theory 

3.1 Relativistic Effects and Vacuum Polarization 

We have been unable to resolve the discrepancy of theory and experiment observed in the recent mea- 
surement [10] . Therefore, we now turn our attention to the evaluation of a few numerically tiny, but still 
important corrections to the Lamb shift, which have not yet been addressed in the literature. 

The authors of Ref. [TUHH] use a value of 205.0282 meV for the first-order Uehling correction SE^ to 
the 2P-2S muonic hydrogen Lamb shift, taken with the exact one-body Dirac wave function for both 
the 2P 1 / 2 as well as the 2S , 1 / 2 states. The relativistic correction thus is 

5E^ } = (205.0282 - 205.0074) meV = 0.0208 meV (3.1) 

for the one-body Dirac theory However, the one-body Dirac theory cannot account for the two-body 
reduced-mass corrections. As shown in Ref. [56j . the full two-body treatment has to be based on a 
generalization of Eq. (|2.7j) to a massive photon, integrated over the spectral function describing the 
virtual electron-positron pairs in the vacuum polarization loop [29] . It is imperative to use the full two- 
body treatment for muonic hydrogen because the parameter governing the relativistic corrections captured 
in the one-body Dirac equation (the fine-structure constant a) is much smaller for muonic hydrogen than 
the parameter characterizing the reduced-mass corrections (ro M /m p ) which are captured in the Breit- 
Pauli Hamiltonian. According to Eq. (25) of Ref. [56], the relativistic correction to the one- loop electronic 
vacuum polarization thus reads 

SE^ = 0.0169 meV, (3.2) 

and we prefer to use this value for our final theoretical evaluations. The additional two-body relativistic 
correction beyond the one-body treatment used in Refs. (TUIEJ] thus is 

SE a = 5E^ - SE^ = -0.0039 meV . (3.3) 



3.2 Self-Energy Vacuum Polarization Corrections 

Because the dominant effect to the Lamb shift is due to the one-loop vacuum polarization potential given 
in Eq. (|2.21[) , it would be highly desirable to carry out the calculation of the one- loop self-energy effect 



10 



Section 3: Advancing Theory 



directly for eigenstates of the full Schrodinger-Uehling Hamiltonian given in Eq. (|2.26[) . If the Uehling 
potential is treated in first order, then the Feynman diagrams are given in Fig. [2] 

We thus consider the modification of the muon self energy due to the Uehling potential, which corresponds 
to the replacement V — > V+V vp in the Schrodingcr Hamiltonian. The effect of high-energy virtual photons 
in the self-energy loops given in Fig. [2] can then be expressed in terms of the Dirac Fi form factor acting 
on the vacuum polarization potential V vp . When rewritten in terms of the noncovariant photon energy 
cutoff e, which is a convenient overlapping parameter in Lamb shift calculations |57| . we have 



p(v P ) _ a 



1> 



V 2 (V + V vp ) 



fm 



10 



(3.4) 



where \ip) is the wave function of the bound state in the full binding potential V 
the Schrodingcr-Coulomb eigenstate, we have in leading order 



Vvp. Denoting by 



4a 4 m% . 
s-^ 6. 



to ■ 



(3.5) 



and thus 



E 



(o) 

H 



4a 5 m 3 R 

o 2 — 3" e 

f 1 



In 



10 
~9~ 



(3-6) 



Numerically, we find that the correction SEjj to the high-energy part due to the vacuum polarization can 
be conveniently expressed in terms of a parameter V^i, 



SE„ = E^ p} - E%> = 



(0) 



11 



ir 2 m?, n 3 



V, 



61 



2c 



(3.7) 



By a numerical diagonalization of the Schrodinger-Uehling Hamiltonian on an exponential grid |39j . we 
find for the Vqi coefficient (the first subscript gives the power of a, the second indicates the power of the 
logarithm), 



V 6 i(2S 1/2 ) = 3.09, V ei (2P 3/2 ) 



-0.023. 



(3.8) 



The anomalous magnetic moment of the electron gives rise to a further combined self-energy vacuum 
polarization correction, 



SE M = 



47rm?, 



1 d(V + K P ) 
r dr 



a ■ L 



r 9r (" 7 



7r 2 m?, 71 3 



Mi 



60 • 



(3.9) 



where ip again is the state in the Uchling+ Coulomb potential, and cj> is the unperturbed state. The matrix 
element Mqq is nonvanishing for P states, 



M 60 (25 1/2 ) = 0, M 60 (2P 1/2 ) 



-0.022. 



(3.10) 



The low-energy part of the muon self-energy, in the Coulomb potential, is 



E 



(0) 



37rm?. n 3 



5m In 



a 2 mpt 



In kg 



(3.11) 



It follows from second-order perturbation theory with time-independent field operators and an upper, 
noncovariant cutoff for the photon energy jSTHSU] . We have numerically calculated the low-energy 6E^ pS> 



part for solutions of the full Schrodinger-Uehling Hamiltonian on an exponential grid and find that it 
can be conveniently expressed as 



SE L = 5E 



(vp) 



SE 



(o) 



V 61 In 



4a 6 7 



37r 2 m?, n 3 



L 



oo 



(3.12) 
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T T— — T 





Figure 3: Light-by-light scattering diagrams. 



where Vqi parameterizes the modification of the logarithmic coefficient and equals the corresponding 
correction from the high-energy part, leading to a cancellation of the e parameter. The coefficient L^o 
gives the modification of the Bethe logarithm due to the Uehling potential. Numerically, we find 



L 60 (25 1/2 ) = 11.28, L 6Q (2P 1/2 ) 



-0.0146. 



(3.13) 



In the total correction 5Eh + 5 El, the auxiliary overlapping parameter e cancels, and we have for the 
combined self-energy vacuum polarization correction, 



8E S 



a m r 



a* 



In 



2a 2 m 



10 



AV 61 + AM® 



-AL 



(in 



-0.0025 meV, 



(3.14) 



when evaluated for the difference of the 2P and 2S states in muonic hydrogen. We have used the defini- 
tions, 



AVer = V 6 i(2P 1/2 ) - V 6 i(2S 1/2 ) , 
AM 60 = M 60 (2P 1/2 ) - M 60 (2S 1/2 ) 
AL 60 = L 6Q (2P 1 / 2 ) - L m {2S 1 / 2 ) . 



(3.15a) 
(3.15b) 
(3.15c) 



Electronic vacuum polarization insertions into the bound muon line have previously been evaluated in 
Eqs. (40) and (45) of Ref. [12] in the leading logarithmic accuracy as 

(3.16) 
(3.17) 



SE [ s °l = -0.006 meV. 



The modification due to our complete treatment thus is 

SE b 



5E svp - 6El° g p = +0.0035 meV. 



3.3 Virtual Light by-Light Scattering 

Finally, we also include the results for the light-by-light scattering graphs shown in Fig. [31 as reported in 
Ref. [SO] j which add up to 

6E LL = -0.00089(2) meV (3.18) 

for the 2P-2S Lamb shift. These entries replace the results for the virtual Delbriick scattering and the 
Wichmann-Kroll term in [TU] (entries # 9 and # 10 in Table I of the supplementary material included 
with Ref. [j"0]). which otherwise add up to 

SE WK+VB = +0.00032(135) meV (3.19) 

The shift is 

SE C = SE Lh - SEwk+vb = -0.0012 meV , (3.20) 
and the corresponding uncertainty of ±0.00135 meV is eliminated from the theoretical prediction. 
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3.4 Zemach Moment Correction 

For S states bound to an infinitely heavy nucleus, the finite nuclear-size (NS) effect is given by 

(i- <»*) 

where (r 2 ) is the mean-square charge radius of the proton, and the a 5 correction to the nuclear size effect 
involves the so-called third Zemach moment (^ 3 )(2) (see Rcf. HO]), which results from the elastic part 
of the two-photon exchange diagrams of bound particle and nucleus, when the both Coulomb photon 
interactions are corrected for the finite-size effect. 

The coordinate-space representation of the third Zemach moment reads 

(r 3 ) (2) = J d 3 r J dV r 3 f E (f- 7) f E (f) , (3.22) 

where Je is the normalized electric charge distribution of the proton, with Jd 3 r f%(r) = 1. 

The third Zemach moment can be obtained, in a model-independent way, from experimentally obtained 
scattering data for electron-proton scattering, in terms of the measured Ge Sachs form factor as 



„3\ 



48 r°° dq 
r-)(2) = — — 



G 2 E {q 2 )-l + \q 2 <r 2 ) 



(3.23) 



where the two subtraction terms ensure the convergence of the integral for small q. From a model- 
independent, careful analysis of the world scattering data performed in Ref. |14) . Friar and Sick obtain 
the value 

(r 3 ) (2) =2.71(13) fm 3 . (3.24) 
We note, in passing, that this analysis excludes the value recently assumed in Ref. [62], which reads 

(r 3 ) (2) = 36.6(7.3) fm 3 . (3.25) 

The magnitude of the estimate put forward in Ref. [62] has independently been called into question [63] . 
For muonic hydrogen, the mass ratio of the orbiting to the nuclear particle roughly is which is still 
small compared to unity. The relative uncertainty of (r 3 ) 2 therefore gives a realistic indication of the 
uncertainty that should be ascribed to the r p term in the formula for the Lamb shift in muonic hydrogen. 

However, as stressed in Ref. [19], the third Zemach moment should not be used directly for an evaluation 
of the nuclear size effect in muonic hydrogen. In Eq. (58) of Ref. |19j . it is shown how to express the two- 
photon exchange diagram, for a finite ratio of m^/rrip, in terms of the proton form factors. Furthermore, 
in Eq. (59) of Ref. [19] , it is shown that in leading order of m^/mp, the third Zemach moment correction 
is recovered. An evaluation of the complete two-photon exchange diagram is then performed for a simple 
dipole model of the nuclear charge distribution, with Ge(—P 2 ) = A 4 /((p 2 + A 2 ) 2 and G M (— p 2 ) = 
\ 9p Ge(—P 2 )- If one assumes the dipole model, then one may relate the third Zemach moment correction 
to the leading finite nuclear-size effect and express (?~ 3 )(2) a s being proportional to r~, where r p is the 
root-mean-square radius of the proton. However, the relation of (?" 3 )(2) and r p is model-dependent, as 
emphasized in a particularly clear manner in Ref. |64j and illustrated on the Gaussian, uniform and 
exponential models for the nuclear charge distribution. 

In order to obtain a realistic estimate for the (elastic part of the) two-photon exchange diagram, which 
captures the model dependence, we thus proceed as follows. We first observe that for the dipole model 
assumed in Ref. [19], the full two-photon exchange diagram yields a correction of +0.018 meV to the 2P- 
2S Lamb shift in muonic hydrogen. With a third Zemach moment based on the dipole model employed in 
Ref. [12], which roughly amounts to (f 3 )(2) = 2.30 fm 3 one obtains +0.021 meV in first order in m M /m p 
[see Eq. (64) of Ref. [19] and the text following this equation], The actual third Zemach moment is 
somewhat larger [see Eq. (|3.24[) ]. but as shown in Ref. [TI5], the reduced-mass effect slightly reduces the 
nuclear-size correction of order a 5 . The ratio of the two results (+0.018 meV versus +0.021 meV), which is 
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equal to 6/7, then gives a realistic estimate for the reduced mass correction to the third Zemach moment 
correction in muonic hydrogen, beyond the factor m M /m p explicitly indicated in Eq. (|3.21j) . However, the 
model-independent value of the third Zemach moment should be used for the final evaluation, as given 
in Eq. (|3 . 24[) . and the correction factor 6/7 should be applied in order to approximately account for the 
reduced-mass effects. This gives a result of 

5E Z = +0.0212(12) meV (3.26) 

for the a 5 correction to the 2P-2S Lamb shift in muonic hydrogen, which replaces the model-dependent 
term 

SE Z = + 0.0347 meV (3.27) 
fm 

used in Refs. [MIS]. 



4 Reevaluation of the Proton Radius 

For convenience, we here recall that in Ref. [10] . the theoretical expression already given above in Eq. 



_ / r 2 r 3 \ 

E th = 209.9779(49) - 5.2262-^ + 0.0347-^r meV, (4.1) 
V fm fm I 

has been used in order to analyze the measured energy interval AE defined in Eq. (|2.ip . 

We here add the two-body treatment of the relativistic correction to the electronic vacuum polarization 
(5E a = —0.0039 meV), the complete result for the vacuum polarization correction to the muon self energy 
(5Eb = +0.0026 meV), and the shift due to the light-by-light scattering graphs (5E C = — 0.0012 meV), to 
obtain a total theoretical shift of 

SE = SE a + 5E b + SE C = -0.0025 meV (4.2) 

for the proton radius independent term. As the uncertainty due to the light-by-light scattering graphs is 
eliminated, the theoretical uncertainty of the proton-radius independent term shrinks from ±0.0049 meV 
to ±0.0046 meV, to which we should add the uncertainty of ±0.0012 meV from Eq. (|3.26j) . The model- 
dependent Tp term given in Eq. (|3.27j) . which is valid only for a dipole model of the nuclear charge 
distribution, is replaced by the mo del- independent term given in Eq. p. 26(1 . Thus, our theoretical pre- 
diction is slightly shifted from Eq. (|2.2p and reads 

E t h = ^209.9974(48) - 5.2262^-^j meV . (4.3) 

A comparison with the experimental result given in Eq. (|2.3p then gives a proton radius of 

r v = 0.84169(66) fm. (4.4) 

This value is only marginally shifted from the value given in Ref. [10], which reads r p = 0.84184(67) fm. 
The difference of the new radius given in Eq. (|4.4p and the CODATA recommended value given above in 
Eq. (jl.ip is 5.0 standard deviations and thus statistically highly significant. 



5 Conclusions 

The purpose of this paper has been twofold. First, in Sec. [2] we have rederived all relativistic and quantum 
electrodynamic contributions to the 2S (F=l) 2P^/ 2 {F^2) transition energy which are relevant on the 
order of the reported discrepancy of theory and experiment |10j . which is 0.31 meV. We have emphasized 
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that the relativistic and quantum electrodynamic corrections relevant on the level of the discrepancy 
represent theoretically well-established corrections. The QED theory of atomic bound states is generally 
recognized as a rather highly developed theory, and the theoretical predictions for muonic hydrogen have 
been obtained as a collective effort of QED theorists [T9H27] . The main effects obtained in the cited 
literature references are independently being verified here. The numerically smaller QED corrections 
listed in Tabic II of Rcf. [35] beyond those treated in Sec. [Hare said to be of "higher order" because they 
are parametrically suppressed (either by higher powers in the mass ratio m^/irip or by higher powers 
in the fine-structure constant a). A correction of a conceivable calculational insufficiency in one of the 
higher-order effects could thus only explain the discrepancy if it leads to a surprising enhancement of the 
corresponding correction that compensates its parametric suppression. 

In Sec. [3] we calculate a number of hitherto neglected effects which contribute to the 2P-2S Lamb shift 
in muonic hydrogen. The most important of these probably is a fully nonperturbative treatment of the 
Uchling correction to the bound-muon self-energy, which is performed using a pseudospectrum of states 
obtained from an exponential grid on a lattice (3S}[35]. We also refer to Ref. [50] for the results on the 
light-by-light scattering graphs, and we present an updated estimate for the a 5 correction to the nuclear 
size effect, which uses a model-independent value for the third Zemach moment [14] . 

We have already stressed that the theory of muonic hydrogen is given, on the level of the discrepancy, by 
only few, simple relativistic and QED corrections not exceeding the two-loop level. In some sense, this 
makes the task easier for theoreticians as compared to the muon anomalous magnetic moment, where 
on the level of the observed 3Aa discrepancy of theory and experiment, a multitude of physical effects 
(three- loop and beyond) contribute (see Refs. [65H681 ) . The observed discrepancy in muonic hydrogen 
should thus be interpreted as a large discrepancy, both in terms of standard deviations as well as in terms 
of its absolute magnitude in energy units. 

The discrepancy is all the more surprising because the spectroscopy of muonic bound states is an estab- 
lished tool for the determination of nuclear radii [5JE] ■ Nuclear radii determined from electron scattering 
and from muonic transitions have been observed to agree on the percent level already in a number of 
experiments, one of the first was reported in 1974 (Rcf. [69| ) . Furthermore, the vacuum polarization 
contributions to the Lamb shift (including two-loop effects) have been verified in other muonic transi- 
tions [70] ■ Possibilities for a clarification of the discrepancy from a theoretical side seem to be somewhat 
limited and will be explored in a following investigation [71] . to the extent possible. 
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